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Combined QCD analysis of e + e data at y/s = 14 to 172 GeV 

O. Biebel a 

a III. Physikalisches Institut A, RWTH Aachen — Physikzentrum, 
Sommerfeldstr., D-53056 Aachen, Germany 

A study of the energy dependence of event shape observables is presented. The strong coupling constant a s 
has been determined from the mean values of six event shape observables. Power corrections, employed for the 
measurement of a s , have been found to approximately account for hadronisation effects. 



1. INTRODUCTION 

Many years of electron-positron annihilation 
experiments have yielded a huge amount of data 
covering a vast region of centre-of-mass ener- 
gies. For example, at the PEP, PETRA, TRIS- 
TAN and LEP colliders annihilation experiments 
have been performed between y/s = 14 and 
172 GeV (see e.g. jjj and references therein). Dis- 
tributions of event shape observables measured 
from hadronic final states clearly show a depen- 
dence on the centre-of-mass energy (e.g. ||). 
This violation of the scaling behaviour is due 
to the running of a s ((i) with the energy scale 

Every measurement of the value of a s from 
event shape distributions requires the treatment 
of impacts of the hadronisation of quarks and glu- 
ons which cannot be calculated perturbatively. 
Usually, phenomenological hadronisation models, 
available from Monte Carlo event generators, are 
employed to correct for hadronisation effects. An 
approach, independent of hadronisation models, 
to account for such non-perturbative effects are 
power corrections H|). In the following, power 
corrections to the mean values of event shape dis- 
tributions are investigated in order to determine 
the strong coupling constant a s from the scaling 
violation that is seen in the energy dependence of 
the mean values. 

2. POWER CORRECTIONS 

Non-perturbative contributions to event shape 
distributions are, on general grounds, expected 



to give rise to corrections proportional to recip- 
rocal powers of y/s. Hadronisation effects, for 
example, should yield a l/\/s correction accord- 
ing to the tube model ||] in which hadrons em- 
anate from a tube spanned in rapidity and trans- 
verse momentum space between quark and anti- 
quark. Another, more formal, theoretical ap- 
proach leading to power corrections deals with 
renormalons g, that are found in the QCD vac- 
uum polarisation at high order perturbative cal- 
culations. Power corrections may also be ex- 
plained in a more phenomenological ansatz j^] 
with the behaviour of the strong coupling at low 
scales GeV. 

The perturbative expression for the energy de- 
pendence of cv s (/i) gives rise to the unphysical 
Landau pole at /x = Aqcd- In general, one would, 
in this energy regime, expect the running such 
that a s (p) remains finite but eventually large for 
all scales down to zero (see e.g. Qj). In the ap- 
proach of H an effective coupling constant 

a p -x(fxi) = ( — a s (fj,)fj, p 
Jo A 4 

is introduced that absorbs all non-perturbative 
details of a s (/i) up to an arbitrary infrared match- 
ing scale fij ~ GeV. Using the effective coupling, 
power corrections to, for example, the mean val- 
ues of event shapes can be calculated. They turn 
out to be simply additive to the perturbative cal- 
culation {T) = (Tpcrt.) + (-T-pow.)- The principal 
structure of the power correction is ^ 

<^pow.) ~ ^"(tQ - lnr (^)' 
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Table 1 

Type of power correction to mean values of event 
shapes. '?' indicates an unknown coefficient. 



observable T 




power terms 


(T) 


-1 


l/V~s 


(C) 


3tt/2 




(M 2 H /s) 


1.0 ±0.5 


1/V~s 


(Bt) 


1.0 ±0.5 


In(Vi)/V5+? • 


(B w ) 


1.0 ±0.5 


In(Vi)/v^+? ' l/>/5 




? 


? • l/s+? • ln(Vi)/s 



• [otp-i(vi) - {a s (v) + ...}] 

in which is a coefficient, p — 1, 2 is the power 
of the term, r = 0, 1 is the power of an 

enhanced power term. All these parameters de- 
pend on the particular event shape. The term in 
curly brackets accounts for the contribution of the 
perturbative expression for the running of a s (/i) 
at n < fii- It is know in next-to-leading order 

The coefficients and parameters of the power 
corrections have been calculated for many event 
shape observables (e.g. ||U), m particular, for 
mean values of thrust T, C parameter, heavy jet 
mass Mh/s, total and wide jet broadening (re- 
spectively, Bt and Bw) and the differential 2-jet 
rate 2/23 using the Durham scheme (cf. Table |l|) 
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Figure 1. Fits to thrust. Top: The solid line is 
the result of a fit of So with a s and fixed. The 
dashed line is the perturbative prediction only. 



Bottom: Fit of a s 
tions. 



and i„ without power correc- 



3. TEST OF POWER CORRECTIONS 

A comparison of the perturbative prediction 
with power corrections with plain perturbative 
predictions has been performed. Figure |] shows 
the results of two fits for the thrust. The first 
fit including the dominating l/\fs power correc- 
tion term obtained x 2 /d.o.f = 1.8 for ao(pi = 
2 GeV) = 0.57 ± 0.01 with fixed values for a s = 
0.122 and 2:^ = 1,0 which is the LEP average from 
event shapes M. The difference to the dashed 
line in Figure pT reveals the size of the power cor- 
rections. The second fit, using the perturbative 
prediction only, resulted in x 2 /d.o.f = 2.4 for 

*If not stated explicitly, a s refers to a s (M z o). 



a s = 0.122 ±0.001 and Xfl = 0.07 ±0.01. The in- 
clusion of power corrections, thererfore, improves 
slighly the description of the data. Thus it allows 
to determine a s without the need for a very small 
renormalisation scale /i = ■ y/s which approxi- 
mates hadronisation effects. 

The two jet broadening measures are expected 
to suffer enhanced power terms. Additionally, 
non-enhanced corrections are expected but the 
coefficients are difficult to predict |^|. Due to this, 
fits have been performed with "enhanced only" 
power corrections, "non-enhanced only" and a 
combination of both. Using fixed values for a s 
and before, the fits obtained x 2 /d.o.f of 

4.6, 3.3, 5.2, respectively, for (Bt) and 11.3, 2.5, 
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Figure 2. Fits to total (top) and wide (bottom) 
jet broadening: The solid lines are the results of 
the fits for ao with a s and x^ fixed. The dashed 
lines are the perturbative prediction only. 



5.5, respectively, for (B\y)- Figure g shows these 
results as solid lines, whereas the dashed line rep- 
resents the plain perturbative prediction. A slight 
preference for "non-enhanced only" power correc- 
tions might be deduced from the fit results for 
(Bw)- However, the data are still not sufficient 
to decide in favour of one of the three variations. 

For (2/23) it is only known that the power cor- 
rection should be quadratic in l/\/s including an 
enhanced term Therefore, the a y coefficient 
has been fitted for, too. Fits have been tried with 
enhanced and non-enhanced power corrections of 
the type l/V* and 1/s. By fixing the values of 



^ 0.04 r 

m 

^0.035 I 
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Fits to mean differential 2-jet rate: 



The solid lines are the results of the fits for 
and cTq with a s and Xu fixed. 



a s and before, all fits obtained x 2 /d.o.f. of 

unity. The a y came out to be of the order 10~ 4 to 
10 -6 and compatible with zero within the errors 
suggesting that power corrections are negligible 
for the mean value of the differential 2-jet rate 
using the Durham scheme. 

4. DETERMINATION OF a s 

The value of a s has been obtained by fitting the 
perturbative calculation plus power corrections to 
the data using a s (x^ ■ M z o) and ao(fii = 2 GeV) 
as parameters of the fit. Since the data still can- 
not yet decide in favour of non-enhanced power 
corrections for the jet broadening, the enhanced 
corrections predicted in Q have been used. For 
the mean value of the differential 2-jet rate using 
the Durham scheme the fit has been performed 
without power corrections as is suggested by the 
results of Section |[ 

Figure |] exemplifies the results of such fits for 
thrust, C parameter and heavy jet mass. For 
comparison the plain perturbative prediction us- 
ing the fitted value of a s is also shown by the 
dashed lines. For the central values of ao and a s , 
x^ = I and Hi = 2 GeV are used. The systematic 
uncertainties for ao are estimated by varying x^ 
between 1/2 and 2 and by varying ajr as stated 
in Table [l]. In addition, the errors of a s include 
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Table 2 

Preliminary results and \ 2 = x 2 /d-0.f. of the fits. 
The errors represent experimental uncertainties, 
x^, and a,f variation. For a s a variation of \xi is 
included. 



T 


a (2 GeV) 


a s (M z o) 


X 2 


(T) 


0.543 ± 0.014 


U.±ZU_ .004 


1.8 


(C) 


0.474 ±0.015 


1 20+ 004 


2.0 


(M%/ S ) 


u - 40 ' -0.077 


n 1 19+0. 006 

U.11Z_ 005 


0.8 


(Bt) 


S42+ 064 

u - J ^-0.036 


nfi+ 007 

U.llD_o.o06 


4.2 


(Bw) 


2fi4+ 048 

U.ZD z i_o.031 


n n 1+O.OO6 
U.lll_o.004 


2.3 


(to) 




n 190+0. 006 

U.1ZO_ .004 


0.9 



a variation of hi by ±1 GeV. The results of the 
fits and the % 2 /d.o.f are given in Table |2|. 

Averaging the six values of a s weighted by their 
total errors, and considering correlations for the 
errors, yields 



a 8 (M z o) = 0.116 ±0.001 



+0.006 
-0.004 



where the first error is the experimental uncer- 
tainty and the second represent the deviations 
due 

quadratically. 
average of a* 



to the variation of Xn, (J-i, and ajr added 
This value agrees with the world 
a -(M z o) = 0.118 ±0.003 |,§. 



5. SUMMARY 

Experiments at e + e _ colliders have measured 
distributions of event shape observables over a 
large centre-of-mass range of yfs — 14-172 GeV. 
The distributions clearly show a scaling violation 
which is expected due to the running of the strong 
coupling constant a s (ii). 

Hadronisation effects contained in these distri- 
butions can be explained by soft gluon phenom- 
ena. Their theoretical treatment yields correc- 
tions proportional to powers of X/s/s in addition 
to the usual perturbative prediction. Power cor- 
rections have been found to be useful approxi- 
mations of hadronisation effects. They have been 
applied to determine a s (M z o) = 0.116io;oo4 irom 
the energy dependence of the mean values of 
thrust, C parameter, heavy jet mass, total and 
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Figure 4. Fits to thrust (top), C parameter 
(middle) and heavy jet mass (bottom): The solid 
lines are the results of the fits for cFo and a s with 
fixed. The dashed line is the perturbative pre- 
diction for the fitted a s only. 
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wide jet broadening, and of the differential 2-jet 
rate using the Durham scheme. The study also 
uncovered that the predicted power corrections 
are not universal for the event shape distributions 
considered. Thus, hadronisation effects cannot be 
explained completely by power behaved correc- 
tions. 
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QUESTIONS 

L. Trentadue, Parma: 

I repeat here, too, the comment I have al- 
ready made to the previous talk by P.A. Movilla 
Fernandez. You have not considered in your anal- 
ysis an event shape variable i.e. the energy-energy 
correlation for which the 1/y/s corrections have 
for the first time been observed. Such a quantity 
has same peculiar characteristics: a) It involves 



experimentally and theoretically a single hadron 
distribution so that hadronisation can be seen as a 
single hadron mechanism instead as a "collective" 
effect; b) hadronisation can be discussed and in- 
cluded analytically without using a Monte Carlo; 
c) a complete next-to-leading resummed expres- 
sion is also known. It would be interesting in my 
opinion to see the results of your analysis for such 
a variable too. 

0. Biebel: 

I know about energy-energy correlations which 
have been investigated already by several experi- 
ments. For this study I focussed on new observ- 
ables for which measurements became available 
by the work of P.A. Movilla Fernandez et al. [Q. 
Additionally, next-to-leading resummed expres- 
sions are known for all of the observables con- 
sidered in this study. I should point out, that I 
used the 0{a 2 s ) perturbative predictions only. 

G. Dissertari, CERN: 

q.q is a universal parameter, so it should be the 
same for all observables. From your analysis it 
looks like this is not the case. Can you comment 
on this? 

0. Biebel: 

Indeed, a universal ao is incompatible with the 
data. I forced ao to be universal in a combined fit 
of thrust, C parameter and heavy jet mass which 
gave rather comparable values for So from the 
individual fits. The combined fit for these three 
yielded a huge % 2 /d.o.f. and also the fit curve 
failed to describe the data. 

G. Dissertori, CERN: 

There is a new prediction for power-law cor- 
rections to the distribution of thrust. Have you 
looked at the distributions, too? 

0. Biebel: 

I have not yet tried this. However, according to 
Yu.L. Dokshitzer and B.R. Webber, the authors 
of this prediction p0| , the calculation is possible 
only for thrust and cannot be transferred to other 
observables. They have fit their prediction in the 
range 0.05 < 1 — T < 0.35 to thrust distributions 
measured at y/s = 14-161 GeV obtaining a sensi- 
ble fit for a s = 0.1185 ± 0.0025 with a x 2 /d-0.f. 
close to unity. 



